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Vibrational nonequilibrium effects in charge transport through single- molecule junctions are in- 
vestigated. Focusing on molecular bridges with multiple electronic states, it is shown that electronic- 
vibrational coupling triggers a variety of vibronic emission and absorption processes, which influence 
the conductance properties and mechanical stability of single-molecule junctions profoundly. Em- 
ploying a master equation and a nonequilibrium Green's function approach, these processes are 
analyzed in detail for a generic model of a molecular junction and for benzenedibutanethiolate 
bound to gold electrodes. 
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Introduction: Charge transport through single- 
molecule junctions, i.e. molecules which are bound 
to metal or semiconductor electrodes, represents an 
interesting and challenging nanoscale nonequilibrium 
problem. Recent experimental advances have allowed 
to study the conductance properties of single-molecule 
junctions and revealed a wealth of intriguing transport 
phenomena [TJ [2j |3j 0]. An important aspect that 
distinguishes nanoscale molecular conductors from 
mesoscopic devices is the influence of the nuclear degrees 
of freedom of the molecular bridge [5j EJ E] • Due to the 
small size of molecules, the charging of the molecular 
bridge is often accompanied by significant changes of 
the nuclear geometry that indicate strong coupling be- 
tween electronic and nuclear (in particular vibrational) 
degrees of freedom [8 . Electronic- vibrational (vibronic) 
coupling manifests itself in vibrational structures in the 
conductance, which have been observed for a variety 
of different systems [5j [6] |9j [TDJ [11] , and may result in 
current-induced vibrational excitation that destabilizes 
the junction and causes local heating p~2j [13]. Fur- 
thermore, conformational changes of the geometry of 
the conducting molecule are possible mechanisms for 
switching behavior and negative differential resistance 
[I4j [15]. The physical mechanisms of many of these 
phenomena, in particular in the resonant transport 
regime, are not well understood. 

In this paper, we investigate vibrational effects in res- 
onant electron transport through single-molecule junc- 
tions with multiple electronic states. The basic mecha- 
nisms of this nonequilibrium problem are illustrated in 
Fig. [I] The interaction of a transmitted electron with 
the vibrational degrees of freedom of the central molecule 
may result in excitation of vibrational quanta (in the fol- 
lowing referred to as emission process, Fig. [l](a)), or de- 
excitationfin the following referred to as absorption pro- 
cess, Fig. [11(b)). For molecular junctions that are domi- 
nated by transport through a single electronic state, these 
basic processes and the resulting conductance properties 
have been studied in detail employing a variety of the- 
oretical methods including master equations, scattering 



theory, nonequilibrium Green's function approaches and 
path integral methods [3[Il[l7j[ia[ll[2ig[^ 

The situation is considerably more complex if multiple 
electronic states of the molecular bridge are involved in 
the transport process. In particular, higher lying elec- 
tronic states can facilitate deexcitation of the vibrations 
via resonant absorption processes associated with the se- 
quential transmission of two electrons (Fig. [I] (c)), thus 
stabilizing the molecular junction. Moreover, vibronic 
coupling results in an effective interaction of electrons 
in different electronic states similar to electron-electron 
coupling in Hubbard-like models. As a result of the in- 
tricate interplay of these different processes vibrational 
and electronic signatures in the conductance of a single 
molecule may become of the same importance. 

Theory: To study vibrationally coupled electron trans- 
port in molecular junctions, we consider a single molecule 
that is covalently bound to two metal leads (L,R) de- 
scribed by the Hamiltonian H=H e \-\-H v [\>. The electronic 
part of the Hamiltonian, 
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FIG. 1: Schematic representation of vibronic transport pro- 
cesses in a molecular junction involving emission (a) and ab- 
sorption (b) of vibrational quanta upon transmission of an 
electron through a single electronic state, as well as emission 
and subsequent absorption associated with sequential electron 
transmission via two different electronic states (c). 
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involves the electronic states of the molecular bridge with 
energies q that are coupled by interaction matrix ele- 
ments Vki to electronic states in the leads with energies 
efe. The vibrational degrees of freedom of the molecular 
bridge are described by 

#vib = ^^4 fl Q + ^2 X ia(a a + a^clci. (2) 

Here, a a is the annihilation operator for a vibration with 
frequency Q a , and X ia denote the corresponding vibronic 
coupling constants. 

Several observables are of interest to analyze the in- 
terplay of the electronic and vibrational degrees of free- 
dom in nonequilibrium transport through a molecular 
junction. Here, we consider specifically the current- 
voltage characteristics J(3>) and the vibrational excita- 
tion (a^ a a a ). To calculate these observables, we em- 
ploy two complementary approaches: a nonequilibrium 
Green's function (NEGF) method [2U [25] as well as 
a master equation (ME) approach. Briefly, both ap- 
proaches are based on the small polaron transformation 
of the Hamiltonian [7], which facilitates a nonpertur- 
bative description of vibronic interaction. The trans- 
formed Hamiltonian H comprises an exactly solvable 
part Hq and a renormalized molecule-lead coupling term 
V = i (ykiXic\,Ci+h.c.) , which involves the shift oper- 
ators Xi = ex.p(^2 a (X ia /Q a )(a a — ah)). ^ is n °ted that 
for models with multiple electronic states Ho includes 
Hubbard-like terms ~ A^A^cJc^ctcj, z^j, which de- 
scribe vibrationally mediated electron-electron coupling. 

In the ME approach, all observables are obtained from 
the reduced density matrix p of the electronic and vibra- 
tional degrees of freedom of the molecule, which is given 
as the stationary limit of the well-established equation of 
motion [TTJ [19] 

Q=-i[H ,p\- dr tr leads {[V, [ V(r) , p p^J ] }-(3) 
Jo 

Here, p^ ads denotes the equilibrium density matrix of 
the leads. In the results reported below, electronic coher- 
ences of the density matrix p are fully taken into account, 
while vibrational coherences are negligible. The current 
is given by I = (2e/h) J °° dr tr{ [F(r), pp£J I}, with 

I = E;, fceL (^4 c ^ " h - c and V ( r ) = e- iH ° T Ve iH ° T . 

The NEGF method is an extension of the approach 
outlined in Refs. 124)1251 to treat multiple electronic states. 
Thereby, the electronic Green's function matrix G is 
determined by the self energy matrix 5] L / R? r/ (t, t') = 
E fce L/R^*^5fc(r,T')<T c X J -(T')X;(T)>, where g k de- 
notes the free Green's function of lead state k and T c 
denotes time-ordering on the Keldysh contour. Vibra- 
tionally induced electron-electron interaction is treated 
non-perturbatively following the strategy of Ref. 26 , The 
calculation of the correlation functions (T c Xj(r f )Xj(r)) 
is based on a cumulant expansion in the dimension- 
less coupling parameters \i a /Q a , which in turn involves 



the electronic Green's function and, therefore, is ob- 
tained in a self-consistent scheme [24] [25]. Within 
the NEGF method, the current is given by I = 
(2e/h) J (dE/(2n)) tr{S<G> - £>G<}. 

In both approaches the molecule-lead coupling is 
treated within (self-consistent) second order perturbation 
theory. The ME approach includes all resonant processes, 
however, it neglects contributions related to co-tunneling, 
and does not account for the broadening of the electronic 
levels induced by the interactions with the leads. These 
processes are taken into account in the NEGF method 
[24] [25], which, as a result, becomes exact for vanishing 
vibronic coupling. On the other hand, the ME approach 
allows an exact treatment of the Hubbard-like terms in 
Ho, while these are approximately accounted for in the 
NEGF method. Thus, both approaches are complemen- 
tary and the good agreement of the results found below 
is a strong indication of their validity. 

Results: First, we consider a generic model system for 
a molecular junction, represented by two electronic levels 
and a single vibrational mode with frequency Q=0.1 eV. 
Both leads are modeled by semi-elliptic bands [18] with a 
band width of 3 eV and an overall molecule-lead coupling 
strength of 0.1 eV. Throughout this article, we consider 
a low temperature of T = 10 K to ensure that vibrational 
excitation is caused solely by nonequilibrium processes. 

Fig. [2] shows the current- volt age characteristic and the 
average vibrational excitation for a moderate vibronic 
coupling, A!/ 2 =0.06 eV, obtained with the ME approach. 
The first electronic level is located at ei=0.15 eV (relative 
to the Fermi energy of the leads) and we consider two 
different locations of the second electronic level, 62- 

As a reference, also the result for a single electronic 
level at e\ is shown. For the latter result, the molecule- 
lead coupling has been enhanced by a factor y/2 to ob- 
tain the same current as in the two-level cases for large 
voltages. Overall, the interaction of the transmitted elec- 
trons with the vibration results in pronounced vibronic 
resonance structures in the current-voltage characteris- 
tic and significant vibrational excitation. As has been 
previously studied in detail for a single electronic level 
[24} [25] , the current exhibits a steplike increase at volt- 
ages $ « 2(ei - Af/fi + raft), ra = 0,1,2,... The cor- 
responding current-induced vibrational excitation also 
shows a step-like monotonous increase with voltage, and 
acquires considerable values. 

The situation becomes significantly more complex if 
another electronic level participates in the transport pro- 
cess. In the current voltage-characteristic, this results in 
a wealth of additional resonance structures. Due to vi- 
brationally induced electron-electron interaction, the res- 
onance, which for a purely electronic transport mecha- 
nism would appear at ~ 2e2, splits into two resonances 
at $ « 2(e 2 -A|/ft) and $ w 2(e 2 - (A| + 2Ai A 2 ) 2 /ft) (in- 
dicated by the dashed lines in Fig.[2|. They correspond to 
the first electronically excited state of the anion and the 
electronic ground state of the di-anion, respectively [27] . 
Vibrational progressions originating from all three elec- 
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FIG. 2: Current and vibrational excitation for a generic 
model system with two electronic levels and a single vibra- 
tional mode obtained for two different locations of the second 
electronic state. The red lines depict results of an approxi- 
mate treatment where the current through the two electronic 
states has been calculated separately and the grey lines give 
results for a single electronic state. 



tronic resonances towards higher voltages can be seen. 
An intriguing effect is the appearance of pronounced vi- 
bronic resonance structures at lower voltages, where the 
electronic level e 2 is still outside the conductance window. 
These structures are related to the absorption processes 
schematically depicted in Fig. [ljc): Upon transmission 
via the lower lying electronic resonance ei, an electron 
may excite vibrational quanta. A subsequently trans- 
mitted electron can absorb these quanta, thus effectively 
lowering the energy of the resonances. This results in 
vibronic structures at voltages $ « 2(e 2 — A^/^ — nfl) 
and $ « 2(e 2 - (A| + 2AiA 2 ) 2 /ft - rafi), n = 0, 1, 2, ... As 
a consequence of the additional vibronic structures, the 
purely electronic resonance steps are barely visible. We 
emphasize that this resonant absorption process is a col- 
lective effect of the two electronic states and cannot be 
described in an approximate treatment, which calculates 
the currents through the two electronic levels separately 
(red lines in Fig. 

The absorption of current-induced vibrational energy 
via higher-lying electronic states has an even more pro- 
nounced effect on the vibrational nonequilibrium distri- 
bution in the stationary state. The results in Fig. [2] 
demonstrate that this process may reduce the vibrational 
energy by more than 50%. Thus, higher lying electronic 
levels stabilize single-molecule junctions with respect to 



FIG. 3: Current and vibrational excitation for a generic 
model system with two electronic states, ei=0.15 eV and 
62=0.6 eV, and different vibronic coupling parameters as in- 
dicated in the legend. 



local heating. 

Another interesting aspect is the influence of the vi- 
bronic coupling strength Ai/ 2 . Fig. 3] shows results ob- 
tained for the two-state model considered above, but dif- 
ferent vibronic coupling parameters. Thereby, the ener- 
gies of the electronic states e 1 / 2 have been adjusted to 
compensate for the different polaron shifts, X^ 2 /Q. In 
the current- volt age characteristic, a larger vibronic cou- 
pling results in more pronounced vibrational step struc- 
tures. The results obtained for the vibrational excitation 
show that, except for regions with resonant absorption, 
the largest current-induced vibrational excitation is ob- 
tained for the system with the weakest vibronic coupling. 
This counterintuitive behavior has been analyzed before 
for the case of a single electronic state [2TJ [22] . In par- 
ticular, it has been shown that the vibrational excitation 
diverges in the limit A^0 if no other vibrational relax- 
ation processes are active. The present results demon- 
strate that this vibrational instability can be considered 
as an 'artifact' of the restriction to a single electronic 
level on the molecule. Vibrational absorption processes 
related to higher- lying electronic states, which will al- 
ways be present in a real molecule, limit the vibrational 
excitation, and thus, prevent this instability. 

The comparison of results obtained with the ME ap- 
proach and NEGF theory in Fig. [3] reveals very good 
agreement for small (data not shown) and moderate vi- 
bronic coupling, and thus, indicates the validity of the 
two approaches. The only major difference is the broad- 
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FIG. 4: Current- volt age characteristic for a BDBT molecular 
junction, with and without coupling to a bath. The dashed 
line depicts results of an approximate treatment [24], where 
the current through the two electronic states has been calcu- 
lated separately. 

ening of the step structures due to tunneling processes. 
For strong vibronic coupling, e.g. \\/2 = 0.12 eV, the 
NEGF approach ceases to be valid due to the cumulant 
expansion employed [24j [25] . 

The results discussed so far were obtained for a generic 
model system of a molecular junction. Finally, we 
consider vibrational nonequilibrium effects in electron 
transport through a realistic molecular junction, where 
a p-benzene-di(butanethiolate) molecule (BDBT) is co- 
valently bound to two gold electrodes. Recent first- 
principle electronic structure calculations for this system 
have shown that the transport properties of this junction 
can be well represented by a model that includes two 
electronic states, localized on the molecular bridge, and 
four vibrational modes [28] . The influence of the remain- 
ing vibrational degrees of freedom of the molecule and 
the phonons of the gold electrodes are represented by a 
coupling to a bath as described in Ref. [24j It is noted 
that, in contrast to the generic model discussed above, 
the two molecular levels that dominate charge transport 
through BDBT are occupied in equilibrium, and thus, 
the prevailing mechanism is hole transport. Correspond- 
ingly, the electronic- vibrational coupling term in Eq. [2] is 

replaced by E a ,iGM(- A ^)(«a + <4) c * c l- 

Fig. [4] shows the respective current- voltage character- 
istics using NEGF theory. Similar as for the model sys- 



tem discussed above, the presence of two vibrationally- 
coupled electronic states results in a splitting of the reso- 
nance corresponding to the lower electronic state into two 
step structures (indicated by the dashed vertical lines). 
The vibronic coupling causes significant current-induced 
vibrational excitation and manifests itself in the current- 
voltage characteristic in a wealth of vibrational structures 
due to emission and absorption processes. As a conse- 
quence, electronic and vibrational features are hardly dis- 
tinguishable. This scenario might be very common for re- 
alistic molecular junctions, which typically involve many 
active vibrational modes. If the coupling to the bath 
is included (solid black line), vibrational relaxation pro- 
cesses compete with current-induced excitation resulting 
in an overall smaller vibrational excitation. As a conse- 
quence, resonant absorption processes are less important 
and electronic resonances in the current-voltage charac- 
teristic are more pronounced. 

Conclusions: Processes due to electronic-vibrational 
coupling can influence the conductance properties and 
mechanical stability of single- molecule junctions pro- 
foundly. The results presented in this paper demon- 
strate that in molecular junctions where multiple elec- 
tronic states participate in the nonequilibrium transport, 
a number of additional vibronic processes have to be con- 
sidered. These include, in particular, resonant absorption 
processes associated with the sequential transmission of 
two electrons as well as vibrationally induced effective 
electron-electron interaction. While the latter results in 
a splitting of electronic resonances, the former may cause 
a wealth of additional vibronic structures in the current- 
voltage characteristic. Moreover, the resonant absorption 
process facilitates vibrational cooling and thus represents 
an important stabilization mechanism of molecular junc- 
tions. Since polyatomic molecules include numerous vi- 
brational modes and often exhibit multiple closely-lying 
electronic states, these findings are expected to be of rel- 
evance for most molecular junctions. 
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